Overcharging by macroions: above all, an entropy effect. 
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Model macroion solutions next to a charged wall show interface true overcharging, charge reversal 
and inversion, and layering. Macroion layering is present, even if the wall or the macroparticle are 
uncharged or if the wall and macroions are like-charged. An effective long-range attractive force 
between the adsorbed macroions is implied. The results are obtained through an integral equation 
theory and a new extended Poisson-Boltzmann theory, and are in accordance with experiments on 
confined macroions and polymer layering. 
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The restricted primitive model (RPM) for an elec- 
trolyte solution includes the two main forces in this sys- 
tem: the long range Coulombic and the short-range re- 
pulsive forces. In RPM the ions are taken to be hard 
spheres of diameter a and charge ezi (e is the protonic 
charge and Zi is the ionic valence) , embedded in a dielec- 
tric medium of dielectric constant e. This model has been 
shown to be in agreement with Monte Carlo (MC) sim- 
ulations and experimental results of bulk and confined 
electrolyte systems Q. 

When a divalent electrolyte, at high concentration, is 
next to a charged wall, the charge of the adsorbed coun- 
terions to the wall overcome that on the wall. This effect 
produces a second layer of ions, where the colons outnum- 
ber the counterions. These effects are known as charge re- 
versal (CR), some times (perhaps improperly) referred to 
as overcharging, and charge inversion (CI), respectively. 
Although these phenomena have been reported, both the- 
oretically j2j and by computer simulations since 1980, 
important implications to protein electrophoresis and 
medicine were later recognized. On the other hand, 
the long-range attraction between confined like-charged 
macroparticles |6) and the adsorption of macroions onto 
oppositely charged [Q or like-charged j|] surfaces have re- 
ceived much attention. The understanding of these phe- 
nomena have been recognized as relevant for the colloid 
science and technology the oil industry, and molec- 
ular self-assembly (e.g., DNA encapsulation) and nano- 
structured films (e.g., polyelectrolyte layering) 

Here, we extend the hypernetted chain/mean spherical 
approximation (HNC/MSA) integral equation to be 
applied to model macroion solutions next to a charged 
wall. The HNC/MSA has been proved to be in good 
agreement with Monte Carlo data for the electrical dou- 
ble layer (EDL) of closely related models . Because 
of the larger macroion's size, this theory is expected to 
be even more reliable than for the simple electrolyte case 
1^,^. The macroparticle is taken as a charged, hard 
sphere of diameter am, concentration pM and valence 



zm, whereas the little ions are modeled by the RPM. 
The wall has uniform surface charge density co- The 
wall dielectric constant is chosen to be equal to that of 
the solvent, to avoid image potentials. The ionic dis- 
tribution, as a function of the distance x from the sur- 
face of the wall, gives the structure of the equilibrium 
EDL, and is expressed in terms of the concentration pro- 
files, Pwiix) = pigwi{x). Pu is the bulk concentration, 
of the ionic species i, and gi„i(a;) is the species i re- 
duced concentration profile (RCP). The HNC/MSA in- 
tegral equations for the RCPs are given by gwi{x) = 
exp[-/3Wi(a;)] = exp [ -/? (eZiV'(x) Ji(a;))]. Wt{x) is 
the potential of mean force, i.e., the effective total wall- 
ion interaction potential. Wi(x) has two contributions: 
the electrostatic potential part, given by the mean elec- 
trostatic potential,?/; (^a;j, plus the short range repulsive 
potential part, due to the ionic size, given by Ji{x). Both 
functions are functionals of pwi{x). /3=1/(KbT), where 
Kb is the Boltzmann constant and T is the absolute tem- 
perature. The ion-ion and the macroion-ion direct inter- 
action potentials are given by a hard-sphere potential 
plus the Coulombic potential. In the limiting case of a = 
HNC/MSA reduces to the integral equation form of a 
new extended inhomogeneous Poisson-Boltzmann (PB) 
theory 0J^,^ for point ions plus macroions, next to a 
charged wall. Since macroions are considered at finite 
concentration, this approach is an improvement to the 
classical PB equation for confined macroions, at infinite 
dilution |^,^, where only two macroions are considered: 
i.e., in our theory macroion-macroion correlations are in- 
cluded. For a—Q and pM —Q, we recover the integral 
equation version of the classical Gouy-Chapman (GC) 
theory for point-ions next to a charged wall . A point- 
ion model (PIM) for an electrolyte solution is like the 
RPM, but a=0. 

We have solved HNC/MSA for several values of Zm, 
om, Pm, Co and salt parameters: z+ : Z-, pi, and a. 
We calculated gwi{x), "^{x) and the eflFective charge den- 
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sity, a{x) = — J pei{y)dy [10,13] . The charge profile 
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in the solution is given by pei{x) = Yy eZmPmgmix), 

m— 1 

where we have omitted the sub-index w, for nota- 
tion simphcity. x=a/2 is the distance of closest ap- 
proach, to the wall, of the small ions. Hence, cto — 

oo 

— / pei{y)dy, by the electroneutrality condition, and 

a/2 

X 

cr'(x) = J" P(.i{y)dy is the charge induced by the wall, 

a/2 

on the fluid, between the wall and the distance x to 
the wall. Hence, a{x) { = a' [x) -\- (Jq) is the effective 
or net charge (wall plus fluid) at the distance x away 
from the wall. cr{x) measures overcharging, CR or CI 
at the interface. The effective electrostatic force on an 
ion is ff{x) = -ez,^^ = ^a{x). Therefore, a{x) 
is not only a measure of the over (under)-charging but 
also of the wall-particle effective electrical force. Fi (x) — 
_dy\^^ — ^^^^^^^^ ^^'^ ^'^^ effective, many-body, 
force between the wall and an ion of species i. Hence, 
F,{x) = /f where //(x) = contains the 

non-electrostatic contributions. The larger the ionic size, 
the larger the confinement excluded volume (CEV) and 
the smaller the accessible volume. In our case, ff{x) has 
the ionic excluded volume contributions. Since both ip(x) 
and Ji {x) are functionals of pei (y) and are in a non- linear 
equation, the charge and size correlations are, in general, 
not independent pO| , p^ and, hence, our theory predicts 
that overcharging, if present, is related to both the elec- 
trical charge and the ionic size. It is an elementary sta- 
tistical mechanics result that the smaller the accessible 
volume the smaller the entropy of the system. Thus, the 
larger the ionic size the lower the entropy of the system. 
Charge correlations of like-charged ions have the effect of 
also reducing the system's accessible volume and, thus, 
of reducing its entropy. 

In all our calculations T=298 K, £=78.5 and pM =0.01 
M. In the RPM calculations a=4.25 A. The 4.25 A 
ionic size approximately corresponds to that of a hy- 
drated ion. In Fig. 0, the macroion has negative charge, 
which is opposite to that of the wall and the divalent 
salt ion. The macroion RCP shows a very strong ad- 
sorption to the wall. A second layer of macroions is ad- 
sorbed, with an intermediate layer of divalent positive 
ions, followed by monovalent negative little ions. A layer 
of positive divalent ions, followed by a monovalent neg- 
ative ions layer, also mediates a third layer of adsorbed 
negative macroions. The macroions of the first layer are 
surrounded by counterions. Subsequent layers of, less 
concentrated, macroions are observed, also mediated by 
layers of positive and negative little ions. Considering 
that the bulk macroion concentration is O.OIM, the local 
macroion concentration at the second peak, «0.035M, 
is not negligible. The first peak is «20M. Hence, large 
macroions, next to a highly charged wall, assemble next 
to the wall. In the inset, at a distance of one macroion 



radius, a deep minimum is observed, corresponding to 
a very strong CR. The maximum located around 15a/2 
show a CI. For x lower than one macroion radius, true 
overcharging of the wall is present: That is, wall's di- 
valent colons, are adsorbed to the wall and their charge 
exceeds that of the wall plus the adsorbed negative lit- 
tle ions. This effect has not been reported before and is 
probably present only in macroions solutions. The effec- 
tive wall electrical field, which is proportional to <t{x), is 
positive before the first layer of macroions and then neg- 
ative, before the second layer. Hence, the electrical force 
is first attractive and then repulsive to negative ions. The 
behavior of the total force on an ion of species i, how- 
ever, is implicit in the RCP, i.e., a gi{x) above (below) 1 
implies that Fi{x) is attractive (repulsive). At x — a/2, 
(j{x) is equal to the wall's charge, whereas for x —^ oo, 
a{x) — > 0, as it should be if electroneutrality is satisfied. 
Lower wall charge density or lower macroion's charge, 
size or concentration produce lower adsorption. The 2:1 
electrolyte solution does not show CR. 

In Fig. 1^, the wall and macroion have negative charge. 
Here, the positive divalent little ions are counterions of 
the wall and macroion. A layer of positive ions, followed 
by negative ions is adsorbed to the wall. Then, a strong 
adsorption of macroions is observed. In the inset, a CR is 
followed by a CI. In Fig. ^ the macroions are responsible 
for the first CR, whereas, in Fig. ^, they are responsible 
for the first CI. The 2:1 electrolyte does show CR. Thus, 
if no macroions are present, monovalent counterions do 
not produce CR, whereas divalent counterions do. The 
first peak corresponds to a 0.2 M macroion concentra- 
tion. This implies a long-range effective attraction, be- 
tween the adsorbed macroions. This attraction is lower, 
as the macroions-wall distance increases or for smaller 
macroions. Our results show that the macroion adsorp- 
tion depends on p^ = p+a'^ + p-a^_ -f pMa\j, i.e., it de- 
pends on the CEV. Hence, adsorption of larger macroions 
implies longer range and more strong effective macroion- 
macroion attraction. These findings are suggestive, in 
relation to experimentally reported attraction between 
like-charged macroions, next to a like-charged wall, where 
/ini size confined macroions show the same behavior [6]. 

In Fig. ^, uncharged, large particles, immersed into 
a 1:1 RPM electrolyte are next to a negatively charged 
wall. A slight CR is observed in the inset, even though 
that the counterions are monovalent. EDL studies for 1:1 
electrolytes show that there is no charge or RCP oscilla- 
tions around 1 Hence, in Fig. ||, the oscillations in 
the little ions RCP are due to size correlations related to 
the macroparticle adsorption to the wall. Thus, macroion 
charge is not necessary to have macroparticle adsorption 
and their presence induce ionic oscillations for 1:1 elec- 
trolytes. 

In Fig. ^ the wall is uncharged. Since the macroion has 
a larger size, aj\/=6.5a, we have increased the macroion's 
charge to have the same macroion's surface charge den- 
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sity. The result is a strong adsorption of macroions to the 
wall. The local macroion concentration, at contact with 
the wall is Ri6.63M. In the inset, strong positive and neg- 
ative true overcharging of the wall is seen. If aM=4.5a or 
Zm=-^0 (not shown) the maximum of the second layer 
of macroions decreases from 1.7, in Fig. ^, to 1.2, but 
the RCP's and a{x) curves are qualitatively equal. If 
no macroions are present, <7{x) — 0,\fx > 0. Since the 
wall is uncharged, the strong macroion adsorption and 
wall overcharging is due to the large CEV, imposed by 
the macroions size. Thus less accessible volume, which 
implies less entropy, impose more order at the interface, 
i.e., more adsorption to the wall. 

In Fig. H, we repeat the calculation of Fig. ||but, now, 
the electrolyte species have zero diameter (PIM) and the 
macroion has a smaller diameter, aA/=4.5a. This case 
corresponds to an extension of the PB theory, where only 
two macroions are considered [||. In the inset, before one 
macroion radius, positive overcharging is present. The 
maximum is at x = aM/2. The maximum of nega- 
tive overcharging is at x ~9.6a/2. If no macroions are 
present, a{x) = 0,\/x > 0. Notice the oscillations in the 
point-ion RCP. This oscillatory behavior is in accordance 
with the exact second moment condition (SMC) of Still- 
inger and Lovett iQ. However, it is well known that 
the RCPs obtained from the PB equation for PIM elec- 
trolytes in bulk, next to a charged wall or around two 
macroions (DLVO theory |^,^), do not show oscillations 
and the macroion-macroion interaction is purely repul- 
sive [0,^. In fact, it is a mathematical theorem that the 
PB equation can not predict and attractive force for two 
like-charged macroions, between them or with the wall, 
at infinite dilution [ p^ . The adsorption of macroions 
to the wall implies an effective attraction between them 
and with the wall. The difference of our Fig. || result 
with the classical PB result is the finite concentration of 
macroions, which implies a proper consideration of en- 
tropy. 

In summary, CR, CI and true overcharging of a wall 
depend on electrostatic interactions and on the CEV, 
which depends on pfp = p+a\ + p-a'L + Pmo^i (larger 
implies more adsorption). On one hand: a) Higher 
wall or macroion charge, enhance adsorption; b) Typi- 
cal hydrated monovalent counterions (e.g., Na^, a «4.25 
A) do not show CR (Fig. |l|), whereas typical divalent 
counterions do (Fig. ||). This is due to a greater effi- 
ciency of divalent counterions to store charge fl^ . On 
the other hand: 1) Overcharging can be present for un- 
charged walls or macroions, or for like charged wall and 
macroions, provided p^ is sufficiently large; 2) Point ions 
can never overcharge a surface [|l|,^,| 13 1; 3) In a pure 
electrolyte solution, larger monovalent counterions, such 
as hydrated Li+ (a «7 A), show CR |l|]. Therefore, 
while that CR, CI and true overcharging result from the 
competition between energy vs. entropy, our results show 
that entropy has a much more important role than has 



been recognized in the past, and give some insight on 
the experimentally found long range effective attraction 
of like-charged macroions next to a wall |^ , polymer lay- 
ering 0decher97, and self-assembled complexes ||]. 
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FIG. 1. Inhomogeneous RCP for macroions in a 2:1 RPM 
electrolyte solution, as function of the distance to the wall 
PM = O.OIM, p+ = 0.7M, p- = l.OM, (TO = 0.272C/m^ 
aM ~ 4.5a, Zm ~ —40, z+ — 2, z- = —1. The solid, dash, 
and dot hues are the macroion (M), negative (-) ion, and 
positive (+) ion RCP, respectively. In the inset the solid line 
is the effective charge density profile, (j{x), as a function of 
the distance to the wall, for the macroion solution, whereas 
the dash line is a{x) for a 2:1 RPM electrolyte (p+ = 0.5M, 
p- = l.OM, z+ = 2, z- = —1), when no macroions are 
present. The sketch roughly reprent the distribution of ions 
indicated by their RCP. 
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FIG. 2. As in Fig. ^ but , cto = -0.272C/m^ In the 
inset, the dash line is a-{x) for a 2:1 RPM electrolyte, when 
no macroions are present, p+ — Q.7M , p_ = 1AM, z+ = 2, 
z- = -1. 
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FIG. 3. As in Fig. |l] but the salt is a 1:1 RPM electrolyte 
solution and p+ = l.OM, p_ = l.OAf, ao = -0.272C/m^ 
aM — 6.5a, Zm — 0, z+ — 1, z- — —1. 
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FIG. 4. As in Fig. |l| but p+ = 0.915M, p- 
(To = O.OC/m^, aM = 6.5a, Za/ = —83. 
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FIG. 5. As in Fig. m but the salt is a 2:1 PIM electrolyte 
solution ao — O.OC/ni , Zm = —40. 
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